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EFFECT OF AERODYNAMIC HEATING ON THE FLUTTER OF 
A RECTANGULAR WING AT A MACH NUMBER OF 2 
By Harry L. Runyan and Nan H. Jones 


SUMMARY 


This paper is concerned with the flutter of a solid wing as affected 
by aerodynamic heating, which can cause a large momentary loss in torsional 
stiffness . Both experimental and analytical studies were conducted and 
good correlation between theory and experiment is shown. 

The cantilever wing which was of solid aluminum-alloy construction, 
was tested "cold" at a Mach number of 2 and did not flutter, but was 
caused to flutter when tested in air preheated to 800° F at a Mach number 
of 2. A large transient loss in torsional stiffness due to aerodynamic 
heating resulted in a short period of flutter. Calculations by the use 
of the theory of Budiansky and Mayers (Journal of Aeronautical Sciences, 
December 1956 ) predicted the time at which the minimum stiffness would 
occur which was very close to the time at which the wing fluttered. 

The aerodynamic theory used for the flutter analysis was the second- 
order theory of Van Dyke (NACA Report H 83 ) . The experimental results 
are compared to a flutter calculation which included the computed loss 
in stiffness due to torsional heating. 


INTRODUCTION 


One of the major structural effects of aerodynamic heating on a solid 
wing is to cause a reduction of torsional stiffness. This loss of stiff- 
ness can be attributed to two causes: first, a change in material prop- 

erties which reduces the modulus of rigidity, and second, a transient 
loss due to thermal stresses set up by a nonuniform chordwise temperature 
distribution which can occur in a highly accelerated flight . The reduction 
due to thermal stresses has been studied by Budiansky and Mayers (ref. l) 
and they have shown that very large decreases in torsional stiffness may 
be encountered for aircraft being rapidly accelerated into high-speed 
flight . 
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Torsional stiffness is one of the primary flutter parameters. For 
some very simple cases, it can be shown that the flutter speed is directly 
proportional to the square root of the torsional rigidity. It is, there- 
fore, obvious that the effect of aerodynamic heating on flutter may be 
important and even at times disastrous. The purpose of this paper is to 
present an experimental flutter result on a solid cantilever wing which 
was tested at a Mach number of 2 in air preheated to 800 F and to compare 
this experimental result with a calculation of the flutter speed and of 
the loss in torsional stiffness. The aerodynamic theory used for the 
flutter analysis was the second-order theory of Van Dyke (ref . 2) . 


SYMBOLS 


A chordwise cross-sectional area, sq ft 

a^,a 2 ,a 3 constants used in equation (3) 


b half chord, ft 

Cp specific heat of air, Btu/lb/°F 

c m specific heat of wing material, Btu/lb/°F 


E 

GJ 

h x 

k 

k 

Li,l 2 ,l 3 

M 


modulus of elasticity, lb/ft 2 

O 

torsional stiffness, lb-in. 

heat-transfer coefficient, Btu/(sq ft)(sec)(°F) 
reduced frequency, boo/V 

conductivity of air, Btu/(sec)(sq ft)(°F/ft) 

,Lk nonlinear aerodynamic coefficients 

2 in eq. (l)) 

Mach number 

x 

first area moment about axis of twist, ft 


(defined 



Prandtl number, Cpi-i/k 
dynamic pressure, lb/sq ft 
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R x Reynolds number, pVx/p 

r radial distance from axis of twist, ft 

t(x) wing thickness, ft 

T temperature at time t, °R 

T s stagnation temperature, °R 

T free-stream static temperature, °R 

00 

T aw adiabatic wall temperature, °R 

V velocity, ft/sec 

x Q axis of rotation measured from leading edge based on chord, 

positive rearward 

x, y, z Cartesian coordinates 

a coefficient of thermal expansion, l/°F 


P = fiP 

7 

‘Hr 


M- 

P 


"m 


1 

ratio of specific heat 
recovery factor 

viscosity, lb-sec/sq ft 
air density, slugs/cu ft 

density of wing material, lb/ft ^ 

axial stress in span direction, lb/sq in. 


t time, sec 

"A time parameter 

a) first torsional angular frequency, radians/sec 

CL 
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Subscripts: 

eff effective 

i initial 


DESCRIPTION OF MODEL AND TESTS 


Model 

The model was constructed of aluminum alloy and had a rectangular 
plan form with a chord of 8 inches and span of llj— inches . The wing had 

a solid cross section which tapered from a 65 AOO 3 airfoil section at the 
tip to a 65A004 at the root. The model was swept back 10° as shown in 
figure 1 in order to raise the divergence speed above the maximum oper- 
ating speed of the tunnel. The model was tested backwards, that is with 
the trailing edge of the 65 A series airfoil acting as the leading edge. 
This was done so that the center of gravity would have a rearward loca- 
tion and thus lower the flutter speed so that it would fall within the 
operating limits of the tunnel. The instrumentation on the model con- 
sisted of two sets of strain gages near the root which were used to meas- 
ure the bending and torsional frequency . 

The model properties are given in the following 1 

Aspect ratio of panel 

Elastic-axis location, percent chord 

Center-of-gravity location, percent chord 

First bending frequency, cps 

First torsion frequency, cps 

Second bending frequency, cps 

Nondimens ional radius of gyration (squared), based on 
half chord 

The wing mass per unit length of span varies linearly from 
O .067 slug/ft at the root to 0.0545 slug/ft at the tip. 


table : 


1.468 

62.5 

57-8 

65 

246 

362 

0.22029 


Wind Tunnel 

The 27- by 27-inch test section of the preflight jet of the Langley 
Pilotless Aircraft Research Station at Wallops Island, Va., was used for 
the test. This tunnel is a blowdown type which exhausts directly to the 
atmosphere. The air could be preheated to approximately 800° F at M = 2. 
The test section and model are shown in figure 2. 
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Test Results 

Two tests were made; the first was conducted with "cold 11 air. The 
second test was made with the air preheated to the maximum temperature 
condition. The wing did not flutter for the cold run. For the hot test, 
the wing began to flutter after being exposed to the airstream for 2 sec- 
onds and continued to flutter for more than 2 seconds and then stopped. 
This phenomenon will be explained in a later section. The total test 
time was approximately 10 seconds . 

The test conditions and flutter results are given in the following 
table : 


Test 

Stagnation 

temperature, 

Op 

Test-section 
density, 
slugs/cu ft 

Test- section 
air velocity, 
ft/sec 

Flutter 

frequency, 

cps 

<1> 

Ib/sq ft 

Cold 

325 

0.00287 

2020 


5855.37 

Hot 

800 

.00204 

2600 

108.6 

6895.2 


ANALYSIS 


This section is concerned with a presentation of the method of flutter 
calculations and of the method of calculating the loss in torsional stiff- 
ness due to aerodynamic heating. 


Method of Flutter Calculations 

The flutter calculations were made using the conventional Rayleigh- 
Ritz type of flutter analysis. Three degrees of freedom were used, 
namely, the uncoupled first bending, second bending, and first torsion. 
The usual flutter determinant as given, for example, in reference 5 was 
used. 


However, instead of employing the more conventional linear unsteady 
aerodynamic theory (ref. 4) in the flutter analysis, the second-order 
theory of Van Dyke (ref. 5 ) was used. This theory takes into account 
the nonlinear effects of airfoil shape and thickness. It has been found 
that, for supersonic speeds, the location of the center of pressure is 
highly dependent on the airfoil shape and, since the location of center 
of pressure with respect to, say, the center of gravity may have very 
large effects on the flutter speed, it was decided to use the more exact 
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nonlinear theory. Since the reduced frequency of the test was small 
(k = O.O87), only first-order terms in frequency were included in the 
nonlinear analysis; however, a check was made for one frequency ratio 
which included third-order frequency terms and no appreciable effect was 
found. The nonlinear aerodynamic coefficients as derived from refer- 
ence 4 are as follows: 


l-l = 0 



= 1 /3 2 - 1 
pk \ B 



Z±L/m\ 6 _a_ + (6m 2 - it) a 
k W 8b 2 3 4 k 8b 2 


M2 


&<■' - 2 *°> 


2+_l/M\ _A_ + A_ 

k U/ 8b 2 p 2 k 8b 2 


^ (1 ' “»> 


r+ l/M\4 


k 2 \3/ 3 2 k 2 8b 2 


-A- (1) 


M,. = 


3k 


4/b 2 -1.\ /2B 2 - l\ , 2] 7+ l/M^My 4 M v ^ 4M 2 /^ _ -.nAxq 

" °\^) ° k W 4b5 3^k 4b3 1 8b 2 


where 

_ (? + 1) *E 

2 p2 


Calculation of Loss of Torsional Stiffness 
Due to Aerodynamic Heating 

The basic theory used in calculating loss in torsional stiffness has 
been intuitively derived in reference 1. Basically, the assumption made 
is that an axial stress Oy "follows the fiber so that in a twisted con- 
dition, a component of cTy can act in such a direction as to introduce a 

twisting moment on the wing. The formula for calculating the effect is 
given in reference 1 and may be written as 

GJ e jj*p = GJ-^ + / Cy-r^dA (^) 

J A 

where a v is the axial stress of an element dA which is located r 

distance from the axis of twist, and the integration is performed over 
the chordwise cross section of the wing. Negative values of 0y indi- 
cate compression and positive values indicate tension. For solid wings, 
such as the one tested, which have most of the mass located near the mid- 
chord, the center portion will not heat up as quickly as the edges. The 
cooler center portion tends to restrain the edges from expanding and, 
thus, causes compressive stresses in the edges which can reduce the 
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effective torsional stiffness. The problem then is to compute the values 
of Oy which are caused by nonuniform heating of the wing. 

The stress cr y at a point x of an airfoil due to a change in tem- 
perature is 


a = E ( a ]_ + a 2 x + a 3 z ) " “ ^i) (3) 

where a-|_, a 2 , and are constants to be determined by boundary con- 

ditions, a is the coefficient of thermal expansion, T is the tempera- 
ture at point x at time r, and Ti is the initial temperature. This 

formula is based on the assumption that plane sections remain plane during 
the deformation. Of course, this assumption is not valid at the tip, 
where the stress must reduce to zero. However, Budiansky and Mayers 
(ref. l) have investigated this tip effect for a free-free beam having 
a double-wedge section. They show that for the aspect ratio of the pres- 
ent wing the change in frequency square is only of the order of 3 percent. 
It is thus evident that the neglect of the tip effect will not materially 
affect the results of this paper. 


The conditions needed for determining the constants are that the 
integral of the stress cr y over the cross-sectional area must be zero 

and that the integral of the first moment of the stress about the axis 
of twist must be zero as follows: 



00 


For a doubly symmetrical airfoil like a symmetrical wedge, both 
a 2 and a^ are zero. For an airfoil having symmetry about one axis, 

say the x axis, then a^ = 0. For the case described herein, the 

65A004 airfoil is symmetrical about the x axis but not about the z axis; 
therefore, a^ and a 2 must be calculated. 
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Temperature Calculations 

The temperature distribution was calculated from the following 
formula : 

T - T ± = (T aw - T ± )(l - e 

where 


and 

T aw = T co(l + n r m2 ) 

This formula is based on the assumption of one-dimensional heat flow 
which implies that there is no chordwise heat flow and that there is no 
temperature gradient normal to the wing surface . 

The temperature distribution is a function of the heat-transfer 
coefficient h x . Because of the relatively rough surface of the airfoil 

in the heat test, it is presumed that the flow across the wing was almost 
entirely turbulent. Therefore, the following turbulent heat-transfer 
formula was used: 



1* - (%)°' 8 Xpr l/3 

where x is the distance from the leading edge, R x is the Reynolds 
number based on x, and Npj. is the Prandtl number. 


Application to a Specific Example 

The foregoing analysis for calculating the change in torsional fre- 
quency has been applied to the present wing. Since no closed analytical 
solution is available for the 65 series airfoil, it was necessary to 
perform the integrations by numerical means. This was accomplished by 
dividing the wing cross section into 18 stations, which were l/20 of the 
chord, and 4 additional stations at the leading edge and trailing edge, 
which were l/40 of the chord. 

The heat-transfer coefficient was calculated for each section by 
using equation (6) . The temperature distribution T - Tj_ was computed 
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from equation (5). The temperature distribution across the chord at 
t = 2 seconds for the tip, midspan, and root is shown in figure 3. 

Note the large change in temperature indicated between the leading edge 
and the 0.6-chord position. 

The following values of the various constants were used in this 
calculation: 


V, 

ft/sec 

2,600 


t s , °R 

1,260 


lb-sec/sq ft . . 

7 X 10"7 


T 0 , °R 

530 

p, 

slugs/cu ft . . 

0.00204 


^r 

0.9 

p m> 

lb/cu ft ... 

168 


k, Btu/(sec)(sq ft)(°F/ft) . . 

9.21 x 10-6 

c m, 

M . 

Btu/lb/°F . . . 

0.21 

2 


N Pr 

0.596 


The loss in torsional stiffness was computed by using equation (2) 
at three spanwise stations - the root, the midspan, and the tip. In 
these calculations the variation of the modulus of elasticity E with 
temperature was taken into account . A plot of the torsional stiffness 
is given in figure 4, where the ratio of the effective stiffness at time 
r to the value at t = 0 is plotted against time . Note that the thinnest 
section, the tip, has suffered a greater loss in stiffness than the thicker 
sections . Since the condition of zero stress at the tip was not satisfied, 
the present calculation overestimates the loss in stiffness at the tip; 
however, it is felt that the tip effect will be relatively small and that 
it can be neglected for the present case. 

With the value of the stiffness computed, the torsional frequency 
and modal shapes were computed by using the iteration procedure of ref- 
erence 5. The bending stiffness was also computed by the use of the pro- 
cedure of reference 5j however, the value of the bending stiffness used 
was calculated at each span station by taking into account the variation 
of E with temperature. 


DISCUSSION OF RESULTS 


The results of applying the method of flutter calculations to the 
present configuration and the calculated operational curve due to aero- 
dynamic heating are shown in figure 5* The velocity coefficient V/bo^ 

is plotted against the frequency ratio c The flutter boundary is 

rather flat for most of the range of frequency ratio but turns up rapidly 
as a frequency ratio of unity is approached. The unstable region is above 


CONFIDENTIAL 


10 


CONFIDENTIAL 


NACA RM L58C31 


the flutter curve. The calculated, operational curve is also shown. The 
numbers shown along the curves indicate the time in seconds. At the 
beginning of the test the wing is in an unstressed condition. The value 
of the flutter speed coefficient V/bo^ is 5*05* the frequency ratio 

(X hl jo^ is 0.262 and is plotted at r = 0 in figure 5. As the wing is 

nonunif ormly heated by the air stream, the torsional and bending frequen- 
cies are changed. The torsion frequency is initially reduced due to the 
stresses resulting from the uneven aerodynamic heating and also to the 
reduction in modulus of elasticity. The maximum change in torsional fre- 
quency occurred at 2 seconds and it had a value of 152 cps or a 58-percent 
change in frequency. The bending frequency at 2 seconds was calculated to 
be 62.8 or a 3.4-percent change from the initial frequency. The flutter 
speed coefficient V^bco a is plotted in figure 5 "the various times as 

indicated up to t = 4 seconds. This operational curve intersects the 
flutter region as indicated at about t = 1 second, and the wing remains 
in the unstable flutter region for about 4 additional seconds at which 
time, the wing, even though hotter, is more evenly heated and has regained 
some of its stiffness. In the experiment, the wing started to flutter at 
about 2 seconds and continued fluttering for slightly over 2 more seconds 
before stabilizing as indicated in the figure. Thus the calculations are 
in fairly good agreement with the incidence of flutter . 


CONCLUDING REMARKS 


This paper has been concerned with the effect of transient aero- 
dynamic heating on the flutter of a solid aluminum-alloy wing. A model 
wing which did not flutter at a Mach number of 2 in air preheated to 
300° F was caused to flutter at a Mach number of 2 when the air was pre- 
heated to a stagnation temperature of 800° F . The flutter is explained 
by the loss of torsional stiffness due to the thermal stresses set up 
as a result of the uneven aerodynamic heating. The flutter speed was 
calculated by using a nonlinear aerodynamic theory based on second-order 
theory of Van Dyke (NACA Report 1183) . The loss of stiffness due to 
aerodynamic heating was calculated and the operational line intersected 
the flutter curve . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 19, 1958. 
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Figure 1.- View of model mounted in test, section. L-9l6lO 
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Figure 3 .- Calculated chordwise temperature distribution at t = 2 sec- 
onds, and M = 2.0. 
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Figure 4.- Calculated loss in torsional stiffness against time at 
M = 2.0. (Tip uncorrected for zero stress . ) 
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Figure 5 . - Effect of aerodynamic heating on flutter . 
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